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Abstract — The Radio Environment Map (REM) provides an ef- 
fective approach to Dynamic Spectrum Access (DSA) in Cognitive 
Radio Networks (CRNs). Previous results on REM construction 
show that there exists a tradeoff between the number of measure- 
ments (sensors) and REM accuracy. In this paper, we analyze 
this tradeoff and determine that the REM error is a decreasing 
and convex function of the number of measurements (sensors). 
The concept of geographic entropy is introduced to quantify this 
relationship. And the influence of sensor deployment on REM 
accuracy is examined using information theory techniques. The 
results obtained in this paper are applicable not only for the 
REM, but also for wireless sensor network deployment. 

Index Terms — Geographic Entropy, Spatial Radio Resource, 
Sensor Deployment, Radio Environment Map 

I. Introduction 

Increases in the number of wireless communication systems 
has created a heterogeneous radio environment where multiple 
Radio Access Technologies (RATs) coexist in the same time 
and space. As a result, User Equipment (UE) with cognitive 
capabilities is crucial for flexible radio resource usage. Mitola 
first proposed Cognitive Radio (CR) in 1998 as a context- 
aware radio technology that can be reconfigured to adapt to 
the environment [ 1 J . The Radio Environment Map (REM) has 
been proposed as a database for dynamic spectrum access 
based on UE location and spectrum usage. It contains multi- 
dimensional cognitive information such as geographic features, 
spectral regulations, equipment locations, radio activity logs, 
user policies, and service providers 0. 

To build a REM, sensors (or UE) must be deployed to 
detect the radio environment. The measurement data from the 
sensors is reported to an REM manager. Several approaches 
have been employed for REM construction. In [3], Grimoud 
et al. used an iterative process to obtain the REM based on 
Kriging interpolation to reduce the measurement data required. 
Riihijarvi et al. [4| developed a probabilistic model for the 
REM which exploits the correlation in the measured data to 
reduce the complexity. And Atanasovski et al. Q produced 
an REM prototype using heterogeneous spectrum sensors. 
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The goal of previous work on REM construction was to 
reduce the number of measurements required and improve 
REM accuracy. There is a tradeoff between the number 
of measurements (or number of sensors), and the accuracy. 
Faint et al. Q examined this relationship using computer 
simulation, and showed that increasing the sensor density 
can increase REM accuracy. However, when the sensors are 
sufficiently dense, the improvement is not significant. In this 
paper, we examine this tradeoff theoretically and determine 
that the relationship between the radio parameter error (REM 
accuracy), is p e — (--=), where M is the number of sensors. 
Besides, we obtain a closed form expression for p e as a 
function of M, which is a decreasing and convex function. 
This verifies the simulation results in J5J. Converse to previous 
approaches, we build the REM by considering the coverage of 
all networks, which is inspired by the Cognitive Pilot Channel 
(CPC) technology in [7 1. Our results are not only applicable to 
REM construction, but also to deployment in wireless sensor 
networks (WSN). 

The rest of this paper is organized as follows. Sensor de- 
ployment and its relationship to REM construction is presented 
in Section II. The analysis of this relationship is provided in 
Section III. In Section IV, we examine the tradeoff between 
the number of sensors and REM accuracy. Section V presents 
some numerical results, and finally some concluding remarks 
are given in Section VI. 

II. Sensor Deployment 

The region is divided into small meshes, which are shown 
as small squares in Fig.Q] Sensors are deployed over the entire 
region, and can be network detectors, spectrum sensing entities 
or just UE. Two sensor deployment schemes are considered, 
one-mesh-one-sensor and random sensor deployment. In the 
one-mesh-one-sensor scheme, a sensor is deployed in each 
mesh randomly. Thus the number of sensors is equal to the 
number of meshes. A sensor measurement is considered to be 
the radio environment for the entire mesh. Thus after gath- 
ering all sensor measurements, the REM can be constructed 
(an example is shown in Fig. |7|ai)). With random sensor 
deployment, the sensors are randomly deployed in the region 
without regard for mesh boundaries. In this case, the majority 
of the sensor measurements in a mesh determines the radio 




environment, and these values are used to construct the REM 
for the region (examples are given in Fig. UI^i) and (ci)). 

III. REM Construction 

A. REM parameters 

Define the binary representation of network k at location 

(x,y) as 



R(k,x,y) 



1 if network k is detected at (x, y) 
otherwise 



(1) 



Radio parameter at a location is characterized by the following 
sum of the binary representations for all networks [0 

T 

/(a;,y) = ^i?(fc,x,2/)x2 fc - 1 (2) 
fe=i 

where T is the number of networks. The radio parameter of 
mesh i is then 

P = argmaxpy (3) 

j 

where is the fraction of the area in mesh i with radio 
parameter j, and Ylj=iPij — !• N = 2 T is the number of 
radio parameters. In Fig. Q] there are 8 radio parameters and 
the radio parameter of mesh 2 is 0. 

The radio parameter error (RPE) of mesh i is defined as 



Pe,i = 1 - maxpjj, 

3 

and the RPE of the entire region is defined as 



Pe 



M 

E 



(4) 



(5) 



where M is the number of meshes. The RPE is not a continu- 
ous and smooth function of pij, thus we define the geographic 
entropy (GE) for convenience. The geographic entropy of a 
mesh is defined as the corresponding uncertainty of the radio 
environment in this mesh. In Fig. Q] we are more certain about 
the radio environment in mesh #1 than that in mesh #2, since 
the radio environment in mesh #2 is more composite. Similar 
to the Shannon entropy [11], the geographic entropy of mesh 
i is defined as 



Hi 



N 

E 



PijlogPij, 



(6) 



and the geographic entropy of the entire region is defined as 

M M N 

H = OLiHi = - E E Pi 3 ^°SPij ( 7 ) 
i— 1 i— 1 n— 1 

where a, is the area fraction of mesh i compared to the area 
of the entire region. For a regular mesh division, such as Fig. 
[U at = jj and 

1 M 

B= Ji y £H i . (8) 

i=i 

B. RPE and GE properties 

In this section, we investigate the geographic entropy and 
radio parameter error, and the relationship between them. 

Theorem 1. The geographic entropy of the entire region is 
0(— h=) —> 0, where M is the number of meshes. 

Proof: The meshes with an impure radio environment are 
distributed along the network boundaries (shown as a solid 
curve in Fig. |2). Denote the length of all these boundaries 
as £, the length of a mesh edge as e, and the length of the 
region edge as L. Then we have M — ( j) 2 , and the number 
of meshes with impure radio environment is upper bounded 



K < 



(9) 



e 

This result is obtained by considering the corresponding 
packing problem along the boundary, which is shown as a 
solid curve in Fig. [2] Moving each point on this line in the 
two normal directions a distance \[2e gives the two dotted 
lines. The area between these lines is 2£v / 2~£- All the meshes 
with an impure radio environment are located between these 
dotted lines, so an upper bound on K is 2^\f2e divided by the 
area of a mesh. An upper bound on the geographic entropy is 
then 

1 »r 1 2V2C, Ar 1 2V2£logiV 

H < —KlogN < — MogiV = -= V 6 — , 

~ M Me VM VS 

(10) 



so that O 



' M 



is an upper bound on H . 



Theorem 2. The RPE of the entire region is O y^fj 0> 
where M is the number of meshes. 



Proof: From 1 — p e .i = rnaxp^ > jf, we have that 



Pe,i < 1 



1 

N' 



The RPE of the entire region is then upper bounded by 

1 2\/2£L 



< 



1 



(11) 



(12) 



which gives the required result. ■ 
Theorems Q] and |2] show the scaling of the GE and RPE as 
a function of M. The relationship between these parameters 
is given in the following theorem. 

Theorem 3. The geographic entropy as a function of the radio 
parameter error is upper bounded by 



A 



(13) 



H <H(p e )+ Pe \og\N-l\ =iP( Pe ) 
Proof: According to Fano's inequality, we have 

Hi < ■ IK , Iog|JV- 1| = l/}( Pe ,i) (14) 

where Hi is the entropy of mesh i. Taking the sum gives 

M i M ( ) ( M \ 

i=l i=l \ i=l J 

(15) 

where (a) is due to Jensen's inequality from the concavity of 

M 

ij)(x). Using p e = 5 E p e ,i> me proof is complete. ■ 

i=i 

The following information theoretic lower bound on the 
entropy as a function of the error probability was given by 
Feder and Merhav ifTUI. 



Lemma 1 ( 11101 ). A lower bound on the entropy h as a function 
of the error probability ir is given by h > 0(tt) where 

02(71" - |) + b 2 



< 7T < h 



< 7T < 



= < 



a,(7r 



-) + h 



with cii = i(i + 1) log (^ji) and bi 



< 7T < 



JV^2 < < N-l 
N-l — — N 



(16) 



= log! 



From [ 10 1, 4>{tt) is an monotone increasing and convex 
function of tt (see Fig. 1 in iflOl . where <f>* corresponds to 
(j> in this paper). Based on Lemma [T] we have the following 
theorem. 

Theorem 4. A lower bound on the geographic entropy of the 
entire region as a function of the radio parameter error p e is 
given by 

H><f>{ Pe ). (17) 



Proof: In Lemma [T] let h = Hi and ir = p e> i so that 

Hi > <j>(p e ,i)- (18) 



The geographic entropy is then 



M M 



(19) 



Since 4>(p e ,i) is a convex function of p e ^, L , a lower bound on 
H is given by 

M ( b ) 1 RI 

H ~mT, ^(P«.*) ^ E^-) = tiPe) (20) 

i—l i=l 

where (b) is due to Jensen's inequality. ■ 
Combining Theorems [3] and |4] we have 

4> (Pe) <H<^ ( Pe ) 

v-- 1 (h) < Pe < r 1 (H) 



(21) 



Thus the geographic entropy is related to the radio parameter 
error, and an increase of (a reduction of) the geographic 
entropy may increase (reduce) the radio parameter error. 

The mesh configuration will also affect geographic entropy 
and the RPE, as shown in the following theorem. 

Theorem 5. If any two meshes are fused, the entropy of the 
entire region will not decrease. 

Proof: Assume meshes 1 and 2 are fused. The radio 
parameter distribution of mesh i is Pn,Pi2, ■ ■ ■ ,PiN> an d the 
area of mesh i is Sj. The radio parameter distribution of the 
fused mesh is 

{Pi,P 2 r-- :Pn} = 

sipn + S2P21 S1P12 + S2P22 siPin + S2P2N 



si + s 2 



Sl + s 2 



•) 1 



si + s 2 



(22) 



As the entropy is concave ifTTl . we have 



1 H(p 11:Pl2 , . . . ,Pin) H -7 H(p2l, P 22, ■ ■ ■ ,P2N) 



S1+S2 S1 + S2 

S1P11 + S2P2I S1P12 + S2P22 SiPin + S2P2N 



< H 



si + s 2 



si + s 2 



si + s 2 



so that 

BiHi + s 2 H 2 < (si + s 2 )H(p' 1 ,p' 2 , . . . ,p' N ). (24) 
Before fusion, the entropy of the entire region is 



H = 



s 1 H 1 + s 2 H 2 , S^' 



S S 
where S is the region area. After fusion, this entropy is 



(25) 



M 

tt, {si +s 2 )H(p' 11 p' 2 ,... 1 p' N ) l=3 

H= + -^. (26) 

From (l24l . H < H' , so the entropy of the region is not 
decreased after fusion. ■ 
Duality provides the following theorem. 

Theorem 6. Any mesh division operation will not increase the 
entropy or the radio parameter error of the region. 



Theorem [6] shows that if some meshes are divided into 
smaller meshes (such as the meshes with composite ra- 
dio propagation environment), and the one-mesh-one-sensor 
scheme is adopted, then the geographic entropy as well as the 
RPE can be reduced. We next examine the tradeoff between 
the number of sensors and the REM accuracy. 

IV. REM Construction Tradeoffs 

In section III, we show that the number of meshes (measure- 
ments, sensors) impact the geographic entropy as will as the 
RPE. In this section, we investigate the tradeoff between the 
number of sensors and the REM accuracy with a near precise 
result. 

A. One-mesh-one-sensor 

If each mesh contains one sensor, the number of sensors 
equals the number of meshes. From Theorem [3] we have 



Pe >V'" 1 (H), 



(29) 



which implies that p e is lower bounded by a function of 
entropy If if)~ l (H) ^ 0, then the sensing error 

can never be reduced to 0. To reduce the probability of error 
requires that p e < f3. From Theorem [2] we have 

(30) 

However, because the upper bounds on p e ^ and K in Theorem 
|2] are loose, the bound in ( f30b is also loose. Therefore, we use 
a probability model to obtain near accurate estimates of K 
and p e i as follow 

Theorem 7. The radio parameter error as a function of M is 

' (31) 



where 



Pc 



In 64 



M 



(32) 



12tt -4V2tanh _1 (l - V2)L 

which is a constant determined by the length of the boundaries 
of all networks £ and the length of entire area 's edge L. 

Proof: Fig. [3] illustrates the boundary of network in a 
unit mesh i with an impure radio environment. This boundary 
can be approximated by a line when M is large. We ignore 
the situation where the boundaries of multiple networks cross 
the mesh, as the probability of this occurring is low when M 
is large. The parameters x and 9 determine a line in Fig. [3] 
where x is the distance between vertex A and the boundary, 
and 9 is the angle between this line and horizontal line. Both x 
and 9 are random variables with probability density functions 
(PDFs) 

/e(0) = -,O<0<^ (33) 



and 



fx{x) 



V2L sin (9 



n< x < 1 J: ii 



(34) 



-^9 



^ 

















- The boundary of 
a network 



A' 



► A mesh, unit square 



Fig. 3. A boundary of network cut a mesh. 



The length of the boundary of network in the mesh is 

X <X\ 
X\ < X < X\ 

The radio parameter error is 



x tan 9 + x cot 9 
i 

cos 9 



J 2 

2 



Pe,i 



sin 29 

tan 9 



X < X\ 
X\ < X < X\ 



£2. 

2 



cos9 2 

where x\, x% and X3 are as shown in Fig. [3] with values 



xi = x 3 



sin 9 



X2 



v^sin (9+ ^ -2sin0 



(35) 

(36) 

(37) 
(38) 



with [*] + = max{0, *}. The expected values of and p e ,i 
are shown in d27l i and d28| i at the top of the next page, and 
the corresponding closed form expressions are 

4V2tanrT 1 (l - V2) 



0.7935 



and 



E[ P e 



it + ln(64) 
12tt 



0.1937 



(39) 



(40) 



where tanh _1 (z) is the inverse hyperbolic function defined as 
tanh _1 (z) = I In j^. We determine the value of K, i.e., the 
number of meshes with an impure radio environment, using 
(|39l . If the first K meshes have an impure radio environment, 
then 



i=l 



(41) 



where (a) is due to the weak Law of Large Numbers (LLN), 
and (6) is from the fact that the meshes with an impure radio 
environment cover all the boundaries of networks. The value 
of K can be estimated as 



K 



(42) 



E[Zi] -4 % /2tanh _1 (l - s/2)e 

If the first K meshes have an impure radio environment 
then the RPE of the entire region is 



1 (c) K w\ l 

i=l 



(43) 



V2L ein(e+ j) 



E[ P e 



(x tan 6 + x cot 6) fx {x)dx - 

V2L B in(f5 + 2f) 



cos 6* 



/e(#)^ (27) 



1 a:* 



2 sin 9 cos ■ 



-Jx(x)dx + 



x tan c 



cos ( 



fx(x)dx fe{9)de (28) 



where (c) is due to the LLN. Substituting the value of K from 
(l42l and the value of I?[p e ,i] from (|40]l in (@4]i gives 



1 7r + ln64 



7T£ 



/M 12tt -4v / 2tanh" 1 (l - \/2)L 



(44) 



where £ is the length of all the boundaries of networks, L 
is the length of the edges of the entire region, and N is the 
number of radio parameters. This is a near precise estimate. 

■ 

Note that this confirms the result p e = O I -4= ) . Similar to 
the derivation of (130b , the number of sensors can be obtained 
using the more precise estimate of p e in d44l . From the bound 
Pe < P, we have 



M > 



(7r + ln64)^ 



M 2 (45) 



12 (-4V2tanh _1 (l - y/2)) L/3 / 

B. Random sensor deployment 

Randomly deploying sensors over the entire region is more 
realistic than one sensor in each mesh. Suppose there are J 
sensors and the region is divided into M meshes. With a 
uniform deployment, the probability that a sensor falls into 
mesh i is -4 Vi Thus the probability that there are no sensors 
in mesh i is 

1 



P0={1 'M 



If J 



kM, then lim po 

M—HX 

that have no sensors is 



J 

(46) 

1V1 J 

k . Thus, the number of meshes 



Mp Q 



M 

,-_,h 



(47) 



The radio parameters for meshes without a sensor are ran- 
domly chosen, so the maximum error probability for an empty 
mesh is still 1 — . An upper bound on the radio parameter 
error is then given by 



1 M 

-Y 

M ^ 



< Pe + ±M P0 (l-i 



(48) 



where p e is obtained from (l44l . Since p* < f3, we have 

(Tr + b.64)^ 



M > 



12(-4v / 2tanh" 1 (l - V2))L(j3 - e - fc (l - ^)) 



A 



(49) 
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Fig. 4. The number of sensors vs. the RPE with five networks. 
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Fig. 5. The number of sensors vs. the RPE with three networks. 



Note that M3 > A/2, i.e., with random deployment the num- 
ber of sensors required to achieve the same error probability 
as the one-mesh-one-sensor scheme is larger. 

Finally, we analyze the situation when a mesh contains 
more than one sensor. In this case, the mesh is divided into 
smaller meshes such that each smaller mesh contains one 
sensor. From Theorem |6] any division operation will not 
reduce the geographic entropy and therefore not reduce the 
error probability. Thus (l48l is still an upper bound on the 
RPE and d49l ) is still a lower bound on the number of meshes 
(sensors) with random deployment. 

V. Numerical Results 

The relationship between the number of meshes (sensors) 
and the radio parameter error (RPE) in shown in Figs. |4]and|5] 
for five and three networks, respectively. This verifies (l44l for 
the one-mesh-one-sensor scheme. The RPE is related to the 
number of sensors M and the length of the network boundaries 
£. As M increases, the RPE decreases, and as £ increases, the 
RPE increases. The value of £ in Fig. |4] is larger than the 
corresponding value in Fig. Thus to achieve the same RPE 
(for example, RPE = 0.04), the number of meshes for five 



- One-mesh-one-sensor scheme: k=1 

- Random deployment scheme: k=3 

- Random deployment scheme: k=4 
Random deployment scheme: k=5 



0.12 0.14 0.16 0.18 

Radio Parameter Error 




(a,) 



Fig. 6. The relationship between the number of sensors and the radio 
parameter error. 
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(a ? , 






networks must be larger than the number with three networks. 
Note that the RPE is a convex function of M, thus when M 
is sufficiently large, the RPE improvement is not significant. 

From d45l l and d49l i, for the same error probability, the 
number of sensors with the one-mesh-one-sensor scheme will 
be smaller that with random sensor deployment. This is 
verified by Fig. [6] As k increases, the number of sensors 
with random deployment approaches the number with the one- 
mesh-one-sensor scheme. This results can be obtained from 
d49l as lim M 3 = M 2 . 

fe— >oo 

Fig. |7] provides three examples of REM construction. In 
Fig. Ilfai), the REM with the one-mesh-one-sensor scheme 
has errors along the network boundaries. This is because 
the radio environment of the meshes along the boundaries is 
impure and may contain errors in the measurement results. 
The radio parameter error of each mesh is illustrated in Fig. 
[7t a 2)- Figs. EJbi) and (ci) illustrate the REM with random 
sensor deployment and k = 1 and k = 2, respectively. 
Note that the results in Fig. EJci) are more accurate than in 
(bi). This is expected since d49l ) indicates that the RPE is a 
decreasing function of k. Figs. |3b2) and (02) show the radio 
parameter error distribution for the REMs in Figs. 0bi) and 
(ci). respectively. A mesh without any sensors is shown in 
red. Note that Fig. |7f C2) contains fewer red meshes than Fig. 
|7jb2). This confirms d47l i. which indicates that the number of 
vacant meshes is a decreasing function of k. In Figs. |7Ja2), 
(b2) and (C2), the meshes that are not red or dark blue denote 
meshes that have an incorrect radio parameter. The number 
of such meshes is lower in Figs. |3b2) and (C2) than in Fig. 
|7ja2). This indicates that random sensor deployment results in 
fewer measurement errors, but there are more meshes with no 
sensors when k > 1. 

VI. Conclusion 

In this paper, we have achieved the relationship between 
the number of sensors and the radio environment map (REM) 
accuracy. The concept of geographic entropy is introduced 
to quantify this relationship. And the influence of sensor 
deployment on REM accuracy is examined using information 
theory techniques. The results obtained in this paper are 



Fig. 7. REM construction results and the corresponding RPE values for the 
networks in Fig. |4] (ai ) corresponds to the one-mesh-one-sensor scheme with 
16 X 16 meshes (sensors); (bi) corresponds to random sensor deployment 
with 16 X 16 meshes and k = 1; and (ci) corresponds to random sensor 
deployment scheme with 16 X 16 meshes and k = 2. The RPE values for 
(ai), (bi) and (ci) are given in (a2), (ba) and (02), respectively. 



applicable not only for the REM, but also for wireless sensor 
network deployment. 
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